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AMENDMENTS TO THE CLAIMS 



1-18. (Cancelled) 

19. (Currently amended) A method for modelling, analyz i ng and/or synth e sizing, processing 
a windowed signal representing sound , the method comprising , by a signal processing 
apparatus, computing simultaneously the frequencies and complex amplitudes from the signal 
using a nonlinear least squares method, whereby the computational complexity is reduced by 
taking into account the bandlimited property of the window resulting in band-diagonal system 
matrices for the computation of the amplitudes and frequency optimization step. 

20. (Previously presented) The method according to claim 19 using a stationary 
nonharmonic signal model according to (Eq. (2)): 



which is a model with K stationary components where each component is characterized by its 
complex amplitude A k and frequency cy fc where w„ is the window of claim 1 ; or 

an harmonic signal model according to (Eq. (3)): 



L k=0 p=0 J 

which is a model with S quasi-periodic stationary sound sources with a fundamental frequency 
o) k , each consisting of S k sinusoidal components with frequencies that are integer multiples of 
co k , in which the complex amplitude of the pth component of the Jdh source is denoted A k , P , and 
where w„ is the window of claim 1 . 

21. (Previously presented) The method according to claim 19 using a nonstationary 
nonharmonic model according to (Eq. (4)): 

r k-ip-i 1 
x n = 3?L ri J]5]^, p (-27r^) p exp(-27ria; fc ^) (4) 





(3) 



fc=0 p=0 
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which is a model with K nonstationary sinusoidal components which have independent 
frequencies co k , in which the amplitudes A kp denote the p-lh order of the k-\h sinusoid, and 
where w„ is the window of claim 1 . 



22. (Previously presented) The method according to claim 20, comprising the computation 
of the spectrum as a linear combination of the frequency responses of the window according to 
(Eq. (11)): 

K-l 

X m = ^A k W(m + u k ) (11) 

A;=0 

for the stationary nonharmonic model, 
or(Eq.(12)): 

S-lS fc -l 

for the harmonic model, 

where the fourier transform of a complex signal results in a spectrum X m> where W(m) denotes 
the discrete time fourier transform of w„ and whereby only the main lobes of the responses are 
computed by using look-up tables. 

23. (Previously presented) The method according to claim 21, comprising the computation of 
the spectrum as a linear combination of the frequency responses of the window according to 
(Eq. (13)): 



X m 



71=0 

K-1P- 



J]^4, p (-2fj^) p exp(-2TO^) 



= EE*- 

k=Q p=0 
K-1P-. 



n=0 
IP 

W(u k + m) 



= EE^4^ fe+m) (13) 

fe=0 p=0 

for the nonstationary model, where the fourier transform of a complex signal results in a 
spectrum X mi where W(m) denotes the discrete time fourier transform of w n whereby only the 
main lobes of the responses are computed by using look-up tables. 
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24. (Previously presented) The method according to claim 20, comprising the step of 
computing the stationary complex amplitudes, by solving the equations (Eq. (19)): 



B i,i b 1 ' 2 

B 2,l B 2,2 



(19) 



where 




Bit = 


N-1 

^ w\ cos(27rw fc 2=^a) cos(27rw z ^a) 

JV-l 


B l d = 


^ wl sin(27TW fc ^i) cos(27rwi^p) 

n=0 


Bit = 


JV-l 

]P wl cos (2iru k 2^22) sm(2irui2=p) 

n=0 


B\i = 


N-1 

n=0 


c\ = 

Cf = 


N-1 

y^ a; n u; n cos(27ru;i^^) 

n=0 
JV-l 

y~] x„tu n sin(27ra;j2rp) 

n=0 



using (Eq. (20)): 

B\i = \*^k + <*)) + \*Vfa-ui)) 

^ 2 = -^(n^+^))-^(^K-^)) 



C, 1 = 3? 



l -J^X m W{m + u l ) 



C ? = -z(~Y,X m W(m + u, l )\ (20) 



such that only the elements around the diagonal of B are taken into account, whereby a shifted 
form B is computed containing only D diagonal bands of B according to (Eq. (27)): 
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B 2 \k = B*f +k _ D (27) 

and Eq. (20), whereby the computation of the Eq. (20) requires the computation of the 
frequency response of the window and the square window denoted by W(m) and Y(m) 
respectively, and solving equation given by Eq. (19) directly from B and C in (Eq. (28)): 

A r = SOLVE(B^,C l ) 

A 1 = SOLVEQP*, C 2 ) (28) 

by an adapted gaussian elimination procedure. 

25. (Previously presented) The method according to claim 20, further comprising the step 
of optimizing the frequencies for the stationary nonharmonic model by solving the equation (Eq. 
(34)): 

HAw = h (34) 
using (Eq. (42)): 

Auj[ = wi —oji 



\A k A l Y"{Cj k + uj l ))-^{A k A]Y" 



Hi k = &(A k AiY'\w k + w l ))-$t(A k AlY''(u> k -u> l )) 

(42) 

such that only elements around the diagonal of H are taken into account, whereby a shifted 
form H is computed containing only D diagonal bands according to (Eq. (36)): 

Hik = Hi t i +k _ D (36) 

and Eq. (42), whereby the gradient h is computed from the residual spectrum R m , where 
Rm = x m - x m denotes the spectrum of the residual r n , and from amplitude Ai and frequencies 
coi, and requires the computation of derivative of the frequency response of the window W'(m), 
whereby the first term of H requires the computation of the second derivative of the frequency 
response of the square window denoted Y"(m), whereby the second term of H is computed 
from the residual spectrum R m , amplitude Ai and frequencies coi, and requires the computation 

of the second derivative of the frequency response W"(m), whereby the parameter X\ allows to 
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switch between different optimization methods and the parameter X 2 regularizes the system 
matrix, and computing the optimization step by solving the system of equations directly on 
Hand h according to (Eq. (37)): 

Aw = SOLVE(H, h) (37) 
by an adapted gaussian elimination procedure. 

26. (Currently amended) The method according to claim 20, further comprising the step of 
optimizing the frequencies for the harmonic signal model, by computing the optimization step 
solving (Eq. (48)): 

HAu = h (48) 

using (Eq. (49)): 

Aloi 
hi 



UJi - Ll>i 

-2- 

N 

S-l P"r 



S,-l /N-l \ 
q=l \m=0 / 



+ J2 qr$t(A p>g Ai, r Y"(qu) p + ru l ) 

/Si-m-i \ 
-A, 5 lp -$l Y, E Rmq 2 A P , q W"(qu; p - m) + 5 lp X 2 (49) 

\( 7 =1 m =0 J 



whereby the gradient h is computed from the residual spectrum R m , where Rm = X m -x m 
denotes the spectrum of the residual r„ and W\m), and from amplitude A { and frequencies ©/, 
and requires the computation of derivative of the frequency response of the window W'(m), 
whereby the first term of H requires the computation of the second derivative of the frequency 
response of the square window denoted Y"(m), whereby the second term of H is computed 
from the residual spectrum R m , amplitude A/ and frequencies oo/, and requires the computation 
of the second derivative of the frequency response W"(m), whereby the parameter Xi allows to 
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switch between different optimization methods and the parameter X 2 regularizes the system 
matrix. 



27. (Previously presented) The method according to claim 19, further comprising the step 
of computing the polynomial complex amplitudes by solving the equation (Eq. (55)): 

B l,l B l,2 

B 21 B 2 ' 2 





A 1 




C 1 




A 2 




C 2 



(55) 



using (Eq. (63)): 



2 [dmP+i 



C lP+q 


■"( 


ClP+q 


= -5 


A 1 

A kP+p 


= 4* 


A 2 

A k.P+p 


= < P 



4[^H_ t „,- ( - 1), il^ a|y( '" , l™,. 



such that only the elements around the diagonal of B are taken into account, whereby a shifted 
form B is computed containing only PD diagonal bands of B according to (Eq. (64)): 

rU _ r1>1 _ rU 

IP+qMP+P - D lP +qt lP +q+ kP+p-{D+l)P+l - n W+q,{k+l-D)P+(p+q-P+l) 

B 2 ' 2 lP+ q ,kP+p = -Bip + g : ;p +9+fc p + p_( jD+1 )p + i = 5jp +gi ( fc+ j_£»)p+(p +g _p + i) (64) 

and Eq. (63), whereby the computation is required of the frequency response of the square 



window and its derivatives a™ p 



Y(m) 



whereby the computation is required of the frequency 



response of the window and its derivatives dLr w ( m \ and solving the equation given by Eq. (55) 
directly from B and C by an adapted gaussian elimination procedure. 
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28. (Previously presented) The method according to claim 24 or 27, comprising a 
preprocessing step which comprises: 

sorting the frequencies to obtain a band diagonal matrix D, determining the number of diagonal 
bands D being defined as the largest k - / for which ~P ^ w k - wz < where u) k and W/ denote 
two frequency values and /3 the width of the main lobe of the frequency response of the 
window. 



29. (Previously presented) The method according to claim 19, further comprising the step 
of computing instantaneous frequencies and instantaneous amplitudes according to (Eq. (69)):, 



$ fc (n) = 2mcj k stfi + i arctan ( ^T? Ak <^~ 2 *^ P \ (69) 

whereby the instantaneous frequency can be used as a frequency estimate for the next 
iteration as expressed in (Eq. (73)): 

(H-l) _ (r) ( , , A k,0 A k,l ~ A kfl A k,l f? oN 
W fc - U k - Kn) ~ 2 ^t- 2 



30. (Previously presented) The method according to claim 19, further comprising the step 
of computing damping factor according to (Eq. (78)): 

Pk « H ^^V^f' 1 (78) 
in case that the amplitudes are exponentially damped. 

31. (Previously presented) The method according to claim 19, where a scaled frequency 
response is used for the analysis of a zero padded window, where wM (' m ~ m o) denotes the 
frequency response of the window of length M and W^{m-n 0 ) = W M (%m-m 0 ) the zero 
padded version of this window up to a length N, and ^'("-"o) = §rW M {§rn-m Q ) the j nverse 

M' = I/V' 

transform of the truncated spectrum to a length N' reducing the window length to 1V1 n , 
resulting in a scaled and zero padded version of the window by computing the inverse 
transform of the scaled frequency response yielding (Eq. (1)): 
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N ' N't i \ 1 \xrM f M ' \ r 0 .(n-nJ)(m-mo). . . 

—wif,(n - n 0 ) = — ^ r(-m - m 0 ) exp(27rz — ^ ) (1) 



32. (Previously presented) The method according to claim 19 for accurate pitch estimation, 
wherein the windowed signal is a sound having a pitch and the method further comprises 
accurately estimating said pitch based on the computed frequencies and complex amplitudes. 

33. (Previously presented) The method according to claim 19, wherein the method is 
applied to determine arbitrary sample rate conversion. 

34. (Previously presented) The method according to claim 19, wherein the windowed signal 
is a sound and wherein noise residual, the amplitudes and the frequencies are encoded in a 
bitstream which is stored, broadcasted or transmitted at a sender side of a 
parametric/sinusoidal audio coder, and a receiver decodes the bitstream back to the 
parameters and synthesizes the sound 

35. (Previously presented) The method according to claim 19 for audio effects whereby 
noise r„, the amplitudes A and the frequencies & are manipulated by an effects processor 
yielding r n> A* anc j o*and synthesized with these modified parameters. 

36. (Previously presented) The method according to claim 19 for source separation, 
whereby sinusoidal components originating from the same sound source are grouped and 
synthesized separately. 

37. (Previously presented) The method according to claim 19 for automated annotation and 
transcription whereby the signal is segmented according to the values of the amplitudes and 
frequencies. 
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